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"spin faster," i.e., increasing the velocity of the circulation, decreasing 
the time required for a given mass to complete the cycle. In either case 
he has increased the energy turn-over per unit of time. Whether, in this 
he has been unconsciously fulfilling one of those laws of nature according 
to which certain quantities tend toward a maximum, is a question well de- 
serving of our attention. 

* Papers from the Department of Biometry and Vital Statistics, School of Hygiene 
and Public Health, Johns Hopkins University, No. 44. 

J Lotka, A. J., Physic. Rev., 24, 1912 (235-238); J. Washington Acad. Sci., 2, 1912 
(2, 49, 66); Science Progress, 55, 1920 (406-417); Proc. Am. Acad. Arts Sci., 55, 1920 
(237-153); these Proceedings, Sci. 6, 1910 (410-415). 

*Proc. Amer. Acad., loc cit., p. 142. 

s See for example Picard, Traite d'Analyse; H. Bateman, Differential Equations, 
1918, p. 245. 

4 Spencer, First Principles, Chapter XXII; Winiarskie, "Essai sur la Mecanique 
Sociale," Revue Philosophique, 44, 1900(113). 

6 At the time of reading proof this project is partially realized. A discussion of the 
applicability of the Le Chatelier principle to systems of the general character here 
considered will appear in a forthcoming issue of the Proceedings of the American 
Academy of Arts and Sciences 

6 Lotka, A. J., London, Phil. Mag., Aug., 1911, p. 353. 



A FORMULA FOR THE VISCOSITY OF LIQUIDS 1 
By H. B. Phillips 
Department op Mathematics, Massachusetts Institute op Technology 
Communicated by A. G. Webster, March 22, 1921 

1. In this paper I obtain for the viscosity of a liquid the formula 

nN h n 

v ~ 2M (»-*)'• U) 

where N is the number of molecules in a mol, h is Planck's constant, M is 
the molecular weight of the liquid in the gas phase, v its volume per gram, 
and n an integer. The quantity 5 is the co-volume as used in the equa- 
tion of state of Keyes 1 

RT A_ 

v-s (v-l) 2 
In all the cases to which I have applied the formula, n = 6 and so (1) takes 

the form 

v(v-s) = Wh/M (3) 

It is to be noted that ZN/M is the number of translational degrees of 
freedom of the molecules in the volume v of the liquid. 

2. To prove equation (1), let x, y, z be rectangular coordinates and sup- 
pose the liquid to flow parallel to the x-axis in such a way that, u being 



*~-T=rrf ( 2 ) 
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the velocity at the point {%, y, z), 

bu /bz=l. (4) 

Let the components of velocity of a molecule parallel to the x and z axes be 
u, w respectively. When a molecule of mass m moves in the positive di- 
rection of the z-axis across the xy-pla.ne, the ^-component of momentum 
transferred across that plane is mu. When it crosses in the negative direc- 
tion the transfer is —mu. Under these conditions the viscosity »/ of the 
liquid is defined as the total ^-component of momentum transferred across 
unit area of the xy-plane in unit time. 2 To find the viscosity we there- 
fore find the number of molecules of each velocity crossing unit area and 
multiply by the average momentum transferred by each. 

Let/ (w) dw be the number of molecules per cubic centimeter with com- 
ponents of velocity parallel to the z-axis between w and w + dw. If the 
molecules did not influence each other, the number of these that would 
cross unit area of the ^y-plane in unit time would be 

wf(w)dw. (5) 

Because of the interference of molecules with each other, this must be re- 
placed by 

v 
— zwf (w) dw. (6) 

11 — o 

To see this we note that if the molecules moved independently, the z- 
component of momentum transferred across unit area of the #y-plane in 
both directions in unit time would be the pressure 

RT/v. (7) 

From equation (2), in the real liquid this is replaced by 

RT/(v-S) * (8) 

Now the mass of the molecule and the temperature (average kinetic energy) 
are the same in both cases. Hence the increased momentum (8) com- 
pared with (7), can only be due to an increase in the number of molecules cross- 
ing in unit time. Since the same law of distribution of velocities (Maxwell's 
Law) is assumed in both cases, the ratio in which the numbers are increased 
must be the same v/(v— 5) for all velocities. Thus the number of mole- 
cules with velocity components between w and w + dw crossing unit area 
of the #;y-plane in unit time is given by equation (6). 

The velocity u of a given molecule parallel to the #-axis will not always 
be equal to the average velocity u of the liquid at that point. There will, 
however, be a series of instants (which we may call collisions) when w will 
be equal to u . I assume that the interval between two such instants will 
be a half period (interval between libration limits) 3 in the sense of the 
quantum theory, and that the molecules can be treated as moving with 
constant velocity between consecutive collisions. Suppose then a mole- 
cule reaching the xy-plane has been moving a time, t, since its last collision. 
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From (4) its momentum parallel to the #-axis due to the motion of the 
liquid at that point will be m w t. 

From (6) the total ^-component of momentum transferred in the positive 
direction across unit area of the zy-plane in unit time will then be 

v f °° 

— I m wHf(w) dw. 

V — o J o 

There will be an equal transfer due to molecules moving in the negative 
direction. Hence 

v f °° 
ij = — I 2mw 2 tf(w)dw. (9) 

v — s J o 

The time t can be considered as the interval between collisions at times 4 

to, h. Since the molecule is assumed to move with constant velocity 

between collisions, 

2mw i t=2 ff 1 mw^dt = nh, (10) 

J U 

where, according to Sommerfeld's theory, n is an integer. Also 



•/: 



00 

f(w)dw=N/2M, (11) 



since it is the number of molecules per gram for which w is positive. Com- 
bining (9), (10), and (11), we get (1) which was to be proved. 

3. Owing to lack of data on the equation of state, the only substances 
on which the formula can at present be tested are carbon dioxide, ether, 
and mercury. Using the values of iV and h given by Birge, 6 

3Mi = 3(6.0594)(6.5543)10-*=.011914, 
and so equation (3) takes the form 

7j(*>-8) = .011914/M. (12) 

Values of 5 for carbon dioxide and ether, 6 tables I and II, were supplied 
by Professor Keyes. At low temperatures the measured viscosities and 
those calculated by equation (12) do not differ by more than the experi- 
mental error. In case of carbon dioxide the temperature 30° is too near 
the critical point (f = 31°) for satisfactory use of the equation of state. 
Above 10° the calculated viscosity of ether is too large, the difference in- 
creasing with the temperature. This may be due to the fact that ether is 
a complex of more than one type of molecule. The equation of state was 
determined on the assumption that each liquid molecule is formed by the 
combination of two gas molecules. This may be substantially true at 
10° and not at 100°. 

To obtain values of v — S for mercury, I make use of the fact that mon- 
atomic substances seem to have constant co-volumes 7 8. Since mercury 
is monatomic in the gas phase, I assume that S is constant or nearly con- 
stant in the liquid phase. Also the term 

RT/{v-8) 
in case of mercury is very large (more than 15000 atmospheres). Hence 
at atmospheric pressures we may neglect p and so write (2) in the form 
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V 



V—8 , n 



where c is constant. Differentiating this equation with respect to v and 
then with respect to T, assuming S constant, we get 

/ dv d 2 v\ /dv\* z/dv\3 

whence 



1 — dv 



dv 

v ->- T i> 



3 + V^^S" 



(13) 



According to Callendar and Moss 8 

Vo = 1+1805553 (ioo) 10 " 8 + 12444 (l5o) 2 10_8 + 2539 (£j 10 " 8 ( M > 

is the ratio of the volume of mercury at t° to its volume at 0°. The values 
of o—5 in table III are obtained by calculating the derivatives of v from 
(14) and substituting in (13). Below 100° the difference of the calculated 
viscosities and those measured is not greater than the difference between 
the values obtained by different observers. At higher temperatures the 
calculated viscosities are consistently smaller than those measured. 
The observations, however, differ greatly, owing probably to oxidation of 
the mercury in contact with the air. The result may be to increase the vis- 
cosity by almost any amount. This work should be repeated with mercury 
of the greatest purity in a vacuum or in contact with some inactive gas. 
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log S=. 48959 



TABLE II 

Ether 
.50981 



' *j(zj-S) = . 0001608 



TEMP. 
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table iii 
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av. 
'F. G. Keyes, "A New Equation of Continuity, 
(323-330). 
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THE ANGULAR DIAMETER OF ALPHA BOOTIS BY THE 

INTERFEROMETER 

By F. G. Pease 

Mount Wilson Observatory, Carnegie Institution op Washington 
Communicated by G. E. Hale, May 2, 1921 

Since the measurement of the diameter of Betelgeuse' in December, 
1920, observations for the determination of stellar diameters have been 
continued with the 20-foot interferometer attached to the upper end of 
the 100-inch Hooker reflector. 

Several definite results have been obtained, and others, less definite 
because of poor observing conditions, have been partially checked by 
comparison with stars whose diameters are known to be too small for 
measurement with this instrument. Let the visibility of the zero fringes 
for any particular separation of the mirrors be reduced by seeing which is 
poor as compared with that on a fine night; it is then hopeless to make 
final measures. If, however, the interferometer is turned to a neighbor- 
ing check star and the fringes are seen, and then turned to the star under 
consideration and they do not appear, the observer is justified in saying 
that there is a definite decrease in visibility. 

The data thus far obtained are as follows: For a Tauri (Aldeberan) 
fringes of gradually decreasing visibility have been observed at 13, 14.5 
and 19 feet. Further observations are necessary between the last two 
positions to determine whether the fringes vanish between these points 
or whether a longer beam will be necessary to obtain a measure of the 
diameter. If the fringes vanish around 16-18 feet, as observations under 
poor conditions lead one to suspect they will, then the fringes observed 
at 19 feet must lie beyond the point of disappearance, on the ascending 
branch of the visibility curve where it rises toward the second maximum. 



